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Abstract 



We study the production and the evolution of quark-gluon plasma expected 
to be formed in ultra relativistic heavy-ion collisions, within the color flux- 
tube model. We introduce the gluonic component in the Boltzmann equation, 
which we solve in the phase space which is extended to include the SU(3) color 
degree of freedom. The color degree of freedom is shown to play a decisive 
role in equilibration, and in fixing the temperature of the plasma. We further 
find that the soft partons that we study here contribute substantially to the 
bulk properties. Finally, the model is shown to provide a detailed picture of 
how the quark-gluon plasma evolves and is driven towards a hydrodynamic 
flow, which starts occuring around 1 fm. 
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I. INTRODUCTION 



It is by now well established from lattice studies |l[ that hadronic matter at sufficiently 
high temperatures (~ 200MeV) and pressures undergoes a transition to the so called quark 
gluon phase involving, i) deconfinement and ii) chiral symmetry restoration, not necesarily 
simultaneously. It also appears that the deconfinement transition is of the first order, and 
that the deconfined phase near the critical temperature is still rather far away from showing 
an ideal gas behaviour. In any case, while such a phase surely did exist in the early universe, 
it is interesting that we may expect to produce this phase in ultra relativistic heavy ion 
collisions (URHIC). In fact it is widely believed that RHIC and LHC will succeed in revealing 
this new phase of matter. 

The deconfined phase is not accessible directly in the accelerator experiments. The sig- 
natures are therefore necessarily indirect, and the prominent ones that have been studied 
are : 1) J/ij) suppression 0, 2) electromagnetic probes such as dilepton and direct photon 
production |§|§], and 3) strangeness enhancement ||. However, a proper diagnostics for 
these signatures in the accelerator produced plasma must involve a study of the evolution of 
the system in all its stages, including its production soon after the two nuclei have collided 
with each other, the equilibration, hydrodynamic expansion leading to cooling, and finally 
hadronization. Note that the above stages are not necessarily mutually exclusive; there could 
be an overlap between production and equilibration, as also between cooling and hadroniza- 
tion. Of particular importance to us here are the production and the equilibration regimes; 
indeed it has been pointed out Offl that J/ip suppression which is generally attributed to 
the equilibrium stage can also arise from the pre-equilibrium stage, an observation which is 
true for strangeness enhancement as well. The contribution of the pre-equilibrium stage to 
the dilepton production and the direct photon production would be even more pronounced, 
and a careful study of the above signals ought to shed light on this stage of the plasma. 

In this paper we discuss the production and the evolution of quark-gluon plasma(QGP) in 
URHIC by employing the Flux-tube Model 0,§ which is a generalization of the familiar Lund 
string model widely used for e + e~ and p—p collisions ||. This model allows for a concurrent 
production and evolution of quark pairs and gluons by a background chromoelectric field. 
Assuming the Bjorken scenario |TU| where a plateau will be seen in the central rapidity region 
in relativistic heavy-ion collisions, the Flux-tube model deals with the baryon free plasma in 
the central region; there would initially be a huge deposit of energy, modeled by the creation 
of a chromo field between the two receding nuclei, which decays spontaneously to produce 
partons. Note that the field configuration which contains this energy is necessarily electric 
like if it has to produce partons. Collisions between the partons and their acceleration ( 
which is again due to the background field) dictate the dynamics of their evolution and 
equilibration. We propose to study this dynamics within the classical Boltzmann equation, 
by explicitly incorporating the production in the source term, acceleration by the background 
chromo-field and also writing down a collision term. 

In a recent paper |TTJ we studied the evolution of a non-abelian qq plasma in the context 
of color flux-tube model, where it was shown that the non-abelian features have a major 
role in the evolution of the system in a manner that was not captured in any of the earlier 



studies ||7|,|q,|T2[1 which were purely abelian in their content. It was also argued that it was 



unlikely that the system would equilibrate instantaneously. However, for simplicity, we had 



ignored there the gluonic component and considered the simpler gauge group SU(2). In this 
paper, we remedy both the drawbacks; we consider the gauge group SU(3) as is appropriate 
for a real plasma. We also include the gluons, which have been completely ignored so far. 
We pay attention to the detailed evolution, and its prediction for different bulk properties 
of the plasma in this paper. Signatures such as dilepton, direct photon and strangeness 
production from this model will be calculated and reported elsewhere. 

There is yet another class of models which describe the evolution of quark-gluon plasma 
such as heavy-ion jet interaction generator(HIJING) [HJ and parton cascade model(PCM) 



14| . These models are governed by perturbative QCD (pQCD). The production and the 
evolution of hard and semi-hard partons are studied by a master rate equation in HIJING and 
by a transport equation in PCM. While these may be reliable in the study of hard partons, 
these perturbative approaches are admittedly insufficient to study the dynamics of the soft 
partons ||13|| , in particular, their production. Being perturbative, they do not incorporate 
any of the non-perturbative aspects such as the formation of the strings and their break up, 
which is studied even in pp collisions. The addition of soft partons to hard and semi- hard 
partons changes the bulk properties of plasma, such as temperature and energy density. Xu 
et al. have observed that this addition leads to an enhanced suppression of J/ip, which can 
be understood to be a consequence of increased number density and a lowered temperature of 
QGP. In our study, we employ the Schwinger mechanism |T3| for particle production, which 



is quintessentially non-perturbative. Even if one were to employ a perturbative version, as 
a time dependent electric field would require it may be noted that an initial electric 
field as a classical saddle point owes its existence to non-perturbative processes, viz, the soft 
gluon exchanges that take place between the two nuclei. On the whole, one may expect that 
the pQCD based studies fll3| , |i4f will be useful in the hard regime, and that the flux tube like 
models will be required to study the soft regime which will become increasingly prominent 
as the system expands and more and more secondaries are produced ]T7| . Indeed, it should 



be possible to develop a unified approach to study the hard and the soft components, say a 
la the approach of Eskola and Gyulassy [TE| who have included the minijets in their so called 



chromoviscous hydrodynamics which is again based on the flux tube model. This study will 
be taken up separately. 

The paper is organised as follows. In section II we set up the transport equations for the 
quarks and the gluons in an extended phase space. Section III contains a description of the 
color flux-tube model, and its incorporation into the transport equations. The numerical 
procedure is presented in section IV and the results are discussed and compared with other 
models in section V. We conclude and summarize the main results in section VI. 



II. TRANSPORT EQUATION IN EXTENDED PHASE SPACE 

In non-abelian theory the color charge is a continuously varying function of time. The 
precession of the color charge (Q a ) obeys Wong's equation p 



dr 

which supplements the Lorentz force equation 



<I<r f abc u,Q b A^ (1) 



Jf- = Q a F a ^u v . (2) 

GST 

Here A atJ- is the gauge potential, and f abc is the structure constant of the gauge group. 
In order to include the color charge in the phase-space, we consider an extended one particle 
phase space of dimension d = 6 + (N 2 — 1), (with N — 3). The extended phase space is taken 
to be the direct sum R 6 © G, where G is the(compact) space corresopnding to the given 
gauge group. In short, in addition to the usual 6 dimensional phase space of coordinates and 
momenta we now have another eight coordinates corresponding to the eight color charges in 
SU(3). In this extended phase space a typical transport equation reads as [pOj]: 



_p M 0" + Q a F; u tfd» + f abc Q a Ayd c Q \ f(x,p, Q) = C(x,p, Q) + S(x,p, Q) (3) 

Here f(x,p,Q) is the single particle distribution function in the extended phase space. 
The first term in the lhs of equation (3) corresponds to the usual convective flow, the second 
term is the non-abelian version of the Lorentz force term and the last term corresponds to 
the precession of the charge as described by Wong's equation. 5 and C on the right hand 
side of equation (3) correspond to the source and collision terms respectively (described 
below). Note that we have to write separate equations for quarks, antiquarks and gluons 
since they belong to different representations of SU(3). 

The term gluonic source merits some explanation here. The classical background field 
that we consider here has, in contrast to the Maxwell field, self interaction. We are interested 
in the stability of the gluonic vacuum ( which is the analogue of the radiation in eletrody- 
namics) against the fluctuations in the classical background field. An adaptation of the 
Schwinger mechanism in QED shows that the fluctuations can indeed produce the gluons, 
i.e. the off-shell classical field can spontaneously produce the on-shell radiative gluonic field 
(see equation 9). There is, therefore, no ambiguity or double counting in this process. The 
source term yields asymptotic gluonic states, and further interaction between the gluons is 
treated separately by a collision term. 



III. THE FLUX-TUBE MODEL 

We briefly review the flux tube model as is appropriate to the context here. In this 
model, the two nuclei that undergo a central collision at very high energies are Lorentz- 
contracted as thin plates. When these two Lorentz-contracted nuclei pass through each 
other they acquire a nonzero color charge (< Q >= 0, < Q 2 >^ 0), by a random exchange 
of soft gluons. The nuclei which act as color capacitor plates produce a chromo-electric field 



between them pl|f22]| . The strength of the field which naturally depends on the strength of 
the color charge residing on the plates cannot be fixed from first principles. We can only fix 
that phenomenologically, say by identifying the field energy with the energy in the central 
region as estimated by Bjorken [|l(]]. This strong electric field creates qq and gluon pairs via 
the Schwinger mechanism which enforces the instability of the vacuum in the presence of 
an external field. The partons so produced, collide with each other and also get accelerated 
by the background field. As described in the previous section, these color charges rotate in 
the color space, a feature which is manifest in this model, but largely ignored in the earlier 



A 



studies [171BIO. Thus the production, collision, acceleration and rotation are implemented 
in one single transport equation, i.e. equation (3). 

It is very difficult to solve the transport equation written above, in general. And we 



have a set of three coupled equations here. As was done in reference ]T1J we make a few 
assumptions. First of all, we admit only those potentials which can be brought to a form 
where the only surviving components are A^ a = (A 01 , A 31 ). This choice restricts to be 
"Maxwell" like, also pointing in the "1" direction in the color space. This restriction is not 
arbitrary as it is known |23| that the non-Maxwellian configurations - where the charges, 
the gauge potential and the fields do not lie in the same direction in the group space - do 
not produce particle pairs, in general. Secondly, we require a boost invariant description 
of the physical quantities |l(J. Accordingly, we demand that the distribution functions also 
be functions of only boost invariant quantities. Finally, we work within the Lorentz gauge 
which is implemented elegantly by the choice A^ a = e^ u d v G a {T), fi, v = 0,3, with all the 
other components zero, r = {t 2 — z 2 ) 1 / 2 is the boost invariant proper time. 

Now we fix the magnitude of the vector charge Q a , which corresponds to the first Casimir 
invariant of SU(3). It is simply the value of the coupling constant. There is also another 
Casimir invariant, viz, d abc Q a Q b Q c , which is also conserved as the QGP evolves. There is, 
however, no way of fixing its value and the experiments presumably impose no restriction 
on its allowed values. In fact, the same holds for lattice analyses as well. For this reason, 
we do not take cognizance of this invariant, and we conveniently resolve the SU(3) charges 
in the polar coordinates: Qi = Q nX=i sinOkCosOi for i ^ 8, and Q$ = Q Ill=i sinOk- 

We now fix the collision term. A collision term can be indeed obtained from pQCD. Apart 
from making the equation hopelessly non-linear, this choice would be good only for the hard 
components, which are not of interest to us here. As aptly pointed by Hung and Shuryak 
24| recently, a microscopic description of collisions gets increasingly cumbersome and also 



unnecessary as more and more secondaries are produced. On the other hand it is admittedly 
true that there is no way to reliably obtain a collision term in the non-perturbative regime. 
So we shall employ a relaxation time approach here, where the relaxation time r c will have 
to be fixed phenomenologically, and in all probability, a posteriori. Recently it is proposed 



25| , |26|| that r c can be local and have a (weak) dependence on r. We take r c to be constant 



here, and further refinements may be incorporated after we have a better understanding of 
the transport phenomena both experimentally and theoretically. 
Within the relaxation approach, the collision term is written as 

c = -PWU-f eq ) ^ (4) 

Here f eq is the distribution function, in local equilibrium. Note that the locality can extend 
to the color space as well apart from space-time. Taking it to be an ideal gas, for simplicity, 
we write, 

2 

Jq ' B exp((p»-Q a At la )u fi /T(T))±l K ' 

where the +(-) sign is to be taken for fermions (bosons). Note that we allow a common 
equilibrium temperature for quarks and gluons. Here m m = (cosht], 0, 0, sinhrj) is the flow 
velocity and rj is the space time rapidity given by tanhrj = z/t. With our choice of potentials, 
it follows that A^u^ = 0, so that the equilibrium distribution function can be written as: 



req 
J < 



exp((^)/T(r)) ± 1 



(6) 



Note also that the temperature depends only on the proper time, in accordance with 
the Bjorken picture. Demanding the same of the distribution functions as well, we require 
that the longitudinal boosts be symmetry operations on the single particle distribution. The 
boost-invariant parameters are, apart from the color charges, 



r={1>-2?) l '*,t=(r,-y),p t = (p&-pl) 1 " 



(7) 



where y = tanh~ 1 (pi/po) is the momentum rapidity. 

The above set of invariant variables also serve to express the source terms for qq and 
gluon pairs which will be obtained by the Schwinger mechanism. The expression for qq 
production (obtained from constant electric field) is written as [[TTJ 



gE cos 6*i 



In 



exp 



2np 



gE cos 6*i 



7T 



1/2 



exp(— a£ 



(8) 



and for gluon pair production the corresponding term in SU(3) is given by the relatively 
stronger term [f20| , pTf 



^(r,e,Pi^i) = (3/2)5 g (r,^,^,^; 



(9) 



We put g — 4 throughout our calculation. 

Now consider the third term f abc A^ a Q b -^f(x,p, Q), in the transport equation (3). Due 
to the restriction of the gauge potentials to the form A^ a = (A 01 , A 31 ), there is an additional 
simplification in the transport equation. In order to see that, note that equation (3) has the 
formal solution 



J 

where £'( r/ ) is given by 

f = sinh" 1 



T — T . 



T, 



p t cosh f T c 



(10) 



7~ . OCOS#i f T , ... 

— smh£ + - / dr"E(r") 

t' p t r' J T > 



(11) 



Clearly, f lhc A 1 Q b -^f(x,p,Q) = 0, as S, f eq and £' only depend on Qi. 



With this final simplification, we get a set of three equations, one each for quark, anti- 
quark and gluon respectively. They are, explicitly, 



d /tanh£ gcos6iE(r)\ d 



pi cosh £ J d£ 



+ Ug = f_Z + MllPhiM 

r c r c p t cosh £ 



(12) 



for quarks and gluons, and 



d_ _ /tanh£ ffcosgiE(r) \ d_ 
dr \ r p t cosh £ / d£ 



f q (r,t,Pt,0i) 



j q _f? , (13) 

r c r c pt cosh £ 

for antiquarks. 

In the process where the field and the particles are present the conservation of energy- 
momentum is expressed in the form: 

d^TZt + W = °< ( 14 ) 

where 

TZt = j pVW, + V q + f 9 )dTdn 7 , (15) 

and 

Tf v = diag(E 2 /2, E 2 /2, E 2 /2, -E 2 /2). (16) 

Here dT = d 3 p/(2ir) 3 p = p t dp t d£ / (2tc) 2 , and dVl 7 is the angular integral measure in the 

color space. The factor 2 in the equation (15) is for two flavors of quarks. Now since energy 
and momentum are conserved in each collision, we have: 



/ 



p u CdTdtt 7 = 0. (17) 



d^Tf + gE(r) J dFdn 7 p v d ^ q 3 + + 4 J dTdtt 7 p u S q + 2 J dTdQ 7 p u S g = (18) 



Taking the first moment of the Boltzmann equation, integrating over the color degrees 
of freedom for f q , f q and f g and making use of equations (14) and (17), we obtain from 
equation (3) 

v d{2f q -_2f q + f g ) 

In the above equation the factor 4 in third term arises because we have two separate transport 
equations for quark and antiquark, each coming with two flavors, and the factor 2 in fourth 
term is due to gg pair production, although there is only one transport equation for gluon. 
Putting v — and 3 in (18) we get two equations, which then yield the following equation 
for the decay of the electric field: 

^ " fdp t p 2 J o °d^nhc£de 1 [2f q - 2f q + /,] + {^/Q)a\E{r)f 2 = 0. (19) 

Here a = <(5/2) exp(0.25/a), a = c(g/2) 5 / 2 ^^ and c = gg^. Finally, C(5/2) = 1.342 
is the Riemann zeta function. 

To solve this equation we fix the form of T(r), by demanding that the particle energy 
density differ negligibly from the equilibrium energy density, in each collision. We then 
relate the proper energy density, which is defined by 

e(r) = J dTdn 7 (p^) 2 (2f q + 2f q + f g ), (20) 

to the temperature by its equilibrium value, whence, 

T(r) = (21) 
We solve the equation (19) numerically to see the evolution of quark-gluon plasma. 
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IV. NUMERICAL PROCEDURE 



The numerical procedure is already discussed in earlier works [|Tl|,|T^]. We use a dou- 
ble self-consistent method. The procedure follows the scheme {T(r) tria i, E(r) tria i} — > 
{/,/, E(t)} -> {f,f} -> T(r) -> ... by repeated use of equations (10), (19), (10), (21), 
which is iterated until there is a convergence to the required degree of accuracy. We have in 
mind the LHC energies, and we have taken the initial energy density as eo = 300GeV/ fm 3 . 
Compared to PCM 0] where the initial particle energy density is taken to be around 1300 
Gev/fm 3 at RHIC, our choice of initial field energy density might appear some what low. 
However, our choice is guided by the estimate that the formation time for a qgp is a fraction 



of a fermi, as we explained in ref [11 1. In any case it would not be very appropriate to 
compare the two initial conditions since the "initiality" is only in a limited sense. Indeed, 
it is the energy density at the location of the receding color plates, and being well in the 
fragmentation region at all later times after the collision, what matters is the energy in the 
central region, corresponding to larger and larger values of r as the system evolves. It will 
be seen that in this region the results that our study yields are not unreasonable, if we make 
a judicious choice for the value of r c . We have studied, in this paper, the results for three 
different values of r c . For hydrodynamic and collisionless limits we have choosen r c = 0.001 
fm and 5.0 fm. We have compared the results of these limitng value of r c to a realistic 
intermediate value r c = 0.2/ra. 



V. RESULTS AND DISCUSSIONS 

The solution of the transport equations following the procedure outlined in the previous 
section allows us to determine directly the temporal evolution of the particle and energy 
densities, temperature and also the rate at which the field energy flows into the particle 
sector. These quantities are of intrinsic interest and are, in principle, amenable to exper- 
imental study via dilepton production. We present these results below. Also of interest 
are the broader questions: when the equilibration sets in, at what time the flow becomes 
hydrodynamic - with or without viscous flow |10|,p7|, etc. 

Then there are features peculiar to our model. Since the color plates are receding away 
from each other, energy is continuously pumped into the field, which subsequently decays 
to produce particles. We study the relative rates at which these two dynamical processes 
proceed in URHIC. This process cannot of course proceed indefinitely. This has already 
been observed by Gyulassy and Csernai in their study of the dynamics of fragmentation 



region pq] . As the plates give up their energy, they decelerate. The deceleration sets limits 
on the times up to which the model is valid. Indeed, the assumption of boost invariance 
breaks down as it is strictly exact only if v p i ate = c. It is safe to assume boost invariance 
so long as v p i ate > 0.9c p9| . A simple estimate shows that this condition, for our choice 
of initial energy density, holds up to ~ 5 — 10/m for LHC energies. There would be other 
competing processes in the fragmentation region, and the Bjorken scenario is probably valid 
upto 3 — 4/m. Keeping this in mind, we have restricted ourselves to (proper) times < 1.5 fm. 
Note that the results presented below will be applicable for RHIC energies only if r < 1/m. 

Finally, the formulation allows us to study an inherently non-abelian quantity, viz, 
< cos 2 6 >, where 9 is the angle between the charge and the field in the group space. Note 



that its abelian counterpart = 1. Wherever possible, we have compared our results with 
the earlier SU (2) study and PCM. We also display the distribution function(s) in the color 
space, which none of the other models can yield, be they flux tube based or pQCD based. 
All the results will be shown for the quarks and the gluons separately. We mention here 
that it is the color Flux-tube model within which one obtains the evolution of the mean 
background chromoelectric field, which is absent, either in HIJING or PCM. The evolution 
of such a background field has a greater impact on the acceleration of the partons present 
in the system, and hence on the observed signatures. For example the cc pair, which is 
produced in early collisions is acted by this back ground field along with cc potential, to 
evolve into a physical J/ip. 



A. Comparison with SU(2) results 



In the earlier work (TT| we had argued that the assumption of the so called abelian 
dominance which was made in a large number of transport studies [[?],[|,[12] does not receive 
any justification from a proper study of a non-abelian transport equation. The present study 
reinforces the same idea, as it is only to be expected. In particular, note that a determination 



of quantities such as the vacuum current or the polarization current |T8| will be particularly 
suspect in view of the fact that the effective charge will only be a fraction of the true charge 
(see below). There now arises the question of the dependence of the results on N if one 
employs the gauge group SU(N) in solving Equation (3). Although we are not in a position 
to make a strict comparison between SU(2) and SU(3) in this paper because the SU(2) study 
had ignored the gluonic terms altogether, it is still useful to see what a limited comparison 
can yield. We shall restrict ourselves to the value r c = 0.2/m to contrast very briefly the 
results of SU(2) and SU(3). 

We shall consider the time dependence of the electric field, the energy density and the 
number density, shown in Figs. 1-3. First of all, it may be noted that the SU(3) quantities 
evolve much more rapidly than their SU(2) counterparts. In fact, while the number density 
and the energy density have attained their maximum value around lfm, at which time the 
field also has considerably decayed, there is hardly any activity in the SU(2) case even upto 
1.5 fm. Indeed, at r ~ lfm, the number density for SU(3) is larger by a factor ~ 10, and 
the energy, by a factor ~ 5, even if we consider only the quark sector. Please note that it 
is not that the energy is merely apportioned between quarks and gluons; the existence of a 
second channel has in no way decreased the flow to the quark sector. 

More significant is the fact that it is the SU(3) flow that shows the required trend towards 
a hydrodynamic flow. The SU(2) coumterparts fail to show any such trend all that way upto 
1.5 fm. We may certainly expect a hydrodynamic flow to occur at much later times, but 
in all likelihood it does not seem to happen at any realistic value of r. We shall discuss the 
SU(3) flow in more detail in the next section. On the whole, the plasma is denser for SU(3) 
due to increased volume in phase space. 



o 



B. Discussion of the results 



We now discuss the results for SU(3) in some detail, for three values of r c - 0.001, 0.2 
and 5 fm. The first (last) choice corresponds to the instantaneous hydrodynamic (collision- 
less) case, in the time scale set by the initial energy density. The quark and the gluonic 
contributions will be shown separately. It is clear from Figs. 4 and 5 that the quarks and 
the gluons show the same behaviour regarding the energy density and the number density 
as functions of proper time. The quark contribution is larger, partly because we have con- 



sidered two flavors, although the source term |27| favors a larger rate for gluons in the color 



space. The situation will probably get reversed if we include the more correct perturbative 
source term for quark [jl(J and gluon production f50 |. In that case, the single and three 



gluon production rates are of the same order as that of two gluons. That the two sectors 
will continue to show the same trend may be reliably assumed. The important common 
feature that Figs. 1, 4 and 5 show is the nature of the evolution at times later than 1 fm. 
The curves suggest an approach to the hydrodynamic flow. Note that a hydrodynamic flow 
would imply that n(r),e(r) and T(r) behave like r^r -4 / 3 and r -1 / 3 respectively. We find 
that the corresponding exponents are —0.7, —1.23 and —0.31 respectively. While the tem- 
perature scaling suggests a close approach to the free flow regime, the other two exponents 
show the presence of drag ||32|| , implying that collisions have not completely ceased. It may 
be expected that full hydrodynamic flow will take over around 2fm. 

First of all let us compare the results for quarks and gluons separately, for r c = 0.2/m 
before presenting the results for different values of r c . In fig-4 we have presented the scaled 
particle energy densites. In contrast to HIJING and PCM, we obtain lower values for energy 
and number densities for gluons than that of quarks and antiquarks together (fig-5). For 
a complete study, a source term for hard parton production is also required and this term 
can be obtained from minijet production at these collider energies, following say ref 



The importance of the non-perturbative contribution may be gauged qualitatively by noting 
that the peak number density for quark-antiquark together(for two flavors) in this model is 
~ 60 /fm 3 in contrast to the PCM value of ~ 120 / fm 3 (for three massless flavors). The 
energy density for quark sectors in PCM is also larger by the same factor, implying that the 
average energy per particle is of the same order in both the cases. 

Let us now consider < cos 2 8 > (where 9 is the angle between color charge and the 
chromoelectric field in color space). This quantity is a good bench mark to characterise the 
'non-abelan'ness of the system, and we present our results in fig-6. This quantity is gauge 
invariant and physical. It may be seen that this value saturates at 9 ~ 7r/4, for both quarks 
and gluons. This value was always unity in an abelian plasma, as 9 = in that case. 

This has a direct impact on the equilibration of the plasma. The equilibration is faster 
around 9 ~ 7r/2 where the background field effect is zero, and is slower at 9 ~ 0. This 
is understood as follows. At larger value of 9, say around tt/2, there is no acceleration of 
partons by the background electric field. Only collisions are present at this angle, and hence 
the rate of equilibration is faster. On the otherhand when the angle 9 = 0, the acceleration 
of partons by background field retards the equilibration of the partons. At any intermediate 
values of 9, the rate of equilibration lies in between these two extreme values. So this average 
value, which is purely due to the non-abelian nature of the theory, has a major role in the 
equilibration of the plasma. This is more clearly shown in the distribution functions which 
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carry the color degrees of freedom explicitly(see fig- 7 and fig-8). 

We conclude the discussion at r c = 0.2/m by making an interesting observation. As 
can be seen in fig-6, the value of < cos 2 6 > fluctutates around its mean value for both 
quarks and gluons, with similar fluctuations in other quantities as well - although it is not 
that pronounced in them. While it certainly indicates that the distribution is approaching 
its equilibrium value asymptotically, the question is whether the fluctuations would persist 
even as the system hadronises. Unfortunately, it is difficult to make any definite assertion 
at this stage. Indeed, as pointed earlier, evolution of the system beyond 3 — 4/m requires 
a full treatment beyond the assumption of the Bjorken scenario in the central region. It is 
also not known when the hadronization exactly sets in. However, if one assumes that the 
fluctuations do indeed persist, it may quite well happen that it will manifest as an anisotropy 
in the parton momentum distribution. This has been pointed out in a study of the related 
'color filament at ion' by Mrowczynski [[31]. It is also possible ||34j| , although we do not know 



in what manner, that these fluctuations in the hadronic phase show up as disoriented chiral 
condensates. It is only a more complete and rigorous study that can settle the status of 
these speculations. 

Now we present the data for the other two values of r c , namely r c = 0.001/m and 
r c = 5/m, corresponding to hydrodynamic and collisionless limits respectively. These results 
are then compared with r c = 0.2 fm. Consider first the number density n(r). This quantity 
is seen to be sensitive to the values of r c . As can be seen from fig-9, n(r) is larger in the 
hydrodynamic limit, where the maximum value is around 120 per fm 3 , at r ~ 1.0/ra. As 
t c is increased to 0.2 fm we find a lesser value, the maximum value being around 80/ fm 3 . 
This trend continues for collosionless limit where particle density is still less. The increase 
in number density as r c is decreased is expected in the context of classical theory. This 
is because, the collision time in any classical non-equilibrium theory depends roughly on 
the inverse of the number density, apart from the other factors. The behaviour of number 
densities in fig-9 for different values of r c in our calculation also reflects the above fact. 

In fig- 10 we present the scaled energy densities for different values of r c . Unlike the 
number density, there is no strong dependence on r c . This means that the average energy 
per particle is different for different values of r c . We find that for r c = 0.001/m, the 
maximum average energy per particle is around 1.0 GeV, whereas it is around 2.0 GeV for 
t c = 0.2/m and 4.0 GeV for r c = 5.0/m. Note that these high energy deconfined partons 
can produce secondary partons, such as strange quarks and also can break a fully formed 
J/ijj as analysed by short- distance QCD in reference |33| . 

Let us consider the electric field. It may be seen from fig-11 that the decay of the field is 
very slow at r c = 0.001/m, compared to r c = 0.2 and 5 fm. Only 30 percent of the field has 
decayed in the hydrodynamic limit in comparison to collisionless and intermediate limits, 
where the decay is very large. At the other end in the collisionsless limit, the decay is at 
a slightly slower rate in comparison to r c = 0.2 fm, implying that the maximal conversion 
rate is for r c around 0.2 fm. As mentioned earlier, one peculiarity of the model is that 
the field energy is continuously created due to the recession of the plates even as the field 
itself decays to produce particles. To study this, in fig-12 we present field energy per unit 
transverse function of ordinary time. For r c = 0.2 fm this energy is much less than 

that at r c = 5 fm and 0.001 fm. This demonstrates that the conversion is indeed dominant 
at r c = 0.2 fm. This is more clearly displayed in fig-13, where we have plotted the ratio of 
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particle energy per unit transverse area to field energy per unit transverse area. 

In fig- 14 we present the evolution of temperature at r c = 0.2 fin. The maximum temper- 
ature we obtain in the color flux-tube model is around 300 MeV. Finally, a brief comment on 
the distribution functions for quarks and gluons with explicit color dependence (see fig-7 and 
8). As discussed earlier, it may be seen that the equilibration is fastest for 9 = n/2, where 
there is no background effect. At 6 = it/ 2 the thermal equilibration occurs at a common 
time t ~ 1.0 fm for both quarks and gluons. 



VI. CONCLUSION 

We have studied the production and the equilibration of a genuinely non-Abelian plasma 
with the color degree of freedom incorporated in both the source term and the background 
term in the transport equation. With the realistic gauge group SU(3) that we have considerd 
here, the distribution functions for quarks and gluons get defined in the extended phase space 
of dimension 14. We have added the gluonic component into the color flux-tube model, which 
was so far absent. 

The role of the color degree and that of gluons was found to be substantial, and non- 
trivial. The plasma is denser than the SU(2) counter part but is comparable to the Abelian 
case. However, it is much cooler than the abelian plasma, where the corresponding temper- 
ature is ~ 800MeV. The value of the effective charge is larger than that of its counterpart 
in SU(2), contrary to naive expectations. The peak energy per particle is ~ 2 GeV which 



indeed is the demarcating scale JT3| between soft and hard processes. In short, the non- 
perturbative aspects of the evolution of QGP offer a rich variety of results which will have 
to be combined with the perturbative studies in order to abtain a complete description of 
the production and evolution of quark-gluon plasma. 
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Figure captions 



FIG. 1. Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for r c = .2fm. The solid line refers to SU(3) case, and the dashed line to SU(2) case. 
FIG. 2. The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for r c = .2fm. The solid line refers to SU(3) case, and 
the dashed line to SU(2) case. 

FIG. 3. The particle number density as a function of proper time (in units of fermi), for 
t c = .2fm. The solid line refers to SU(3) case, and the dashed line to SU(2) case. 
FIG. 4. The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for r c = .2/m. The solid line refers to total energy density, 
upper dashed line to quark plus antiquark energy density, and lower dashed line to the gluon 
energy density. 

FIG. 5. The particle number density as a function of proper time (in units of fermi), for 
r c = .2/m. The solid line refers to total number density, upper dashed line to quark plus 
antiquark number density, and lower dashed line to the gluon number density. 
FIG. 6. < cos 2 6 > as a function of proper time (in fermi) at r c = .2/m for quark (solid 
line) and gluon (dashed line) 

FIG. 7. f q / fg q as a function of proper time at p t = 300MeV, £ = and r c = 0.2 fm, for 
different values of 9. Solid line corresponds to 9 = n/2. 

FIG. 8. f g / fg 9 as a function of proper time at p t = 300MeV, £ = and r c = 0.2 fm, for 
different values of 9. Solid line corresponds to 9 = it/ 2. 

FIG. 9. The particle number density as a function of proper time (in units of fermi), for 
r c = .2/m (solid line), for r c = 0.001/m (upper dashed line) and for r c = 5/m (lower dashed 
line). 

FIG. 10. The particle energy density scaled w.r.t the initial field energy density as a 
function of proper time (in units of fermi), for r c = .2/m (solid line), for r c = 0.001/m 
(upper dashed line) and for r c = 5/m (lower dashed line). 

FIG. 11. Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for r c = .2/m (solid line), for r c = .001/m (upper dashed line) and for r c = 5/m (lower 
dashed line). 

FIG. 12. The field energy/unit transverse area(in units of GeV//m 2 ) for r c = 0.2 fm (solid 
line), for r c = 0.001 fm (upper dashed line) and for r c = 5 fm (lower dashed line), as a 
function of ordinary time (in fermi). 

FIG. 13. The ratio of particle energy per unit transverse area to field energy per unit 
transverse area for r c = 0.2 fm (solid line), for r c = 5 fm (upper dashed line) and for t c = 
.001 fm (lower dashed line), as a function of ordinary time (in fermi). 
FIG. 14. Evolution of temperature as a function of proper time (in fermi), for r c = 0.2 fm. 
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